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Abstract. Thick disks are faint and extended stellar components found around several disk galaxies including our 
Milky Way. The Milky Way thick disk, the only one studied in detail, contains mostly old disk stars («10Gyr), 
so that thick disks are likely to trace the early stages of disk evolution. Previous detections of thick disk stellar 
light in external galaxies have been originally made for early- type, edge-on galaxies but detailed 2D thick/thin 
disk decompositions have been reported for only a scant handful of mostly late-type disk galaxies. We present in 
this paper for the first time explicit 3D thick/thin disk decompositions characterising the presence and properties 
(e.g. scalelength and scaleheight) for a sample of eight lenticular galaxies by fitting 3D disk models to the data. 
For six out of the eight galaxies we were able to derive a consistent thin/thick disk model. The mean scaleheight of 
■"■». the thick disk is 3.6 times larger than that of the thin disk. The scalelength of the thick disk is about twice, and its 

central luminosity density between 3-10% of, the thin disk value. Both thin and thick disk are truncated at similar 
radii. This implies that thick disks extend over fewer scalelengths than thin disks, and turning a thin disk into a 
thick one requires therefore vertical but little radial heating. All these structural parameters are similar to thick 
disk parameters for later Hubble-type galaxies previously studied. We discuss our data in respect to present models 
for the origin of thick disks, either as pre- or post-thin-disk structures, providing new observational constraints. 

a 

Key words. Galaxies: photometry - Galaxies: bulges - Galaxies: structure - Galaxies: fundamental parameters - 
Galaxies: evolution - Galaxies: formation - Galaxies: individual: ESO311-012, NGC1596, NGC2310, NGC3564, 
NGC3957, NGC4179, NGC4521, NGC5047 

1. Introduction tion with larger scaleheight compared to the inner, domi- 
nating 'thin disk'. 

The knowledge of the detailed distribution of stars in 

galaxies is of fundamental importance to address the for- There are three distinct families of hypotheses for its 

mation and evolution of those systems. To a first approx- creation. The first group considers the thick disk as a 

imation, a disk galaxy can be described by a set of dis- separate entity produced in an early phase of enhanced 

tinct stellar entities: a disk population, a bulge compo- star formation during the initial proto-galactic collapse 

nent, and a stellar halo. Deep surface photom e try of ex- li^the ELS scenario, jEggen, Lynden-Bell, ft ftffiiW 

ternal early-type galaxies feursteinl ll979albO iTsikoudl U2S2 iGjhnori [1984 i Burkert, Truran, & Hensleri U222D ■ 

HorfllllQScft and later elab orate measurements in our own Another family of models regards the thick disk as 

Cxalaxv llGilmore fc Reidl . Il98i revealed the need for an an extension (by dynamical heating) of the thin disk. 

additional compo nent of stars. This was called 'thick disk' The y assume that after the lmtial colla P se a11 § as set " 

teursteid.ll979dY since it exhibiting a disk-like distribu- tles down mto the S alactic P lane and starts formin S 

stars. On this thin stellar disk a variety of constant 

or viol ent heating mechani s ms could act: spiral densit 



o 



- 1—^ 

x 



or violent Hearing mecnamsms could act: spiral density 
wave s ijBarbanis fc Woltierl Il967t ICarlberg k. SellwoodL 
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* Based on observations obtained at the European Southern 19851). encounters with giant molecular clouds (e.g. 
Observatory, Chile ISnitzer fc Schwarzsch ild. 1953t lLacevMl984h . scattering bv 
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Table 1. Global parameters of the observed lenticular galaxies: (i) Principal na me, (2) right ascension , (3) de clination, 
(4) RC3 coded Hubble-type, and the (5) Hubble parameter T are taken from Ide Vaucouleurs et ail l)l99ll) . The (6) 
diameter in arcminutes, the (7) heliocentric radial velocities, and the B-Band absolute magnitude (10) are taken from 
LEDA. According to the heliocentric radial velocities corrected for the Local Group infall into the Virgo cluster (8) 
from LEDA, we estimated the (9) distances following the Hubble relation with the Hubble constant from the HST key 
project of H = 72 km s^Mpc -1 (|Freedman et alll200lL 



massive black holes l Lacev fc Ostrikerl Il985l). ener gy 



irra 

put by accretion of satellite galaxies llCarnev et allll98£ ; 
Quiim Hernauist. fc Fullagarl Il993t IVelazauez fe Whit" , 
199 91 lAguerri et all l200l[K ~or b ar bend i ng ins tabilities 



Raha et all Il99l|) . For example, iGnedinl l)2003|) recently 
used N-body simulations to show that tidal heating in 
a cluster is sufficient to thicken stellar disks by a fac- 
tor of 2-3. This kinematic heating and vertical expan- 
sion will lead to a significant morphological transforma- 
tion of a normal spiral galaxy into a lenticular. The third 
model suggests that thick disks are mostly made of de- 
bris material from accreted satellites. Recen t cosmologi- 
cal N- body+SPH galaxy formation models of lAbadi et alJ 
( 2003]) locating the thick disk formation before z « 1 find 
that more than half of the thick disk stars are actually 
tidal debris from disrupted satellites. Therefore the thick 
disk is not a former thin disk thickened by a minor merger. 
To decide which of these hypotheses could explain the 
thick disk phenomenon best we need first a more general 
and complete statistic of thick disk properties. Naturally, 
these are rather global ones for external galaxies whereas 
our particular position in the Milky Way makes it possible 
to determine much finer details. 



Since t he work of Tsi k oudi and Burstein I BursteinL 
Il979albllcl iTsikoudil Il979l Il980h it appears well known 
that thick disks are quite common in SO galaxies. 
However, none of the more recent detections, except 



for two galaxies in Ide Griis fc van der KruitJ l|l996j) and 
a short remark : in Ide Griis fc Peletierl |)199t1). quantify- 
ing detailed parameters such as the ratio of thick to 
thin disk scaleheight or scalelength, is actually made 
in SO galaxies. All these galaxies are of later Hubble 
type. In addition, we have not found a detailed 2D 
thin/thick disk decomposition for any SO galaxy in the 
literature. Subsequent numerical decompositions deal- 
ing with SO galaxies after the pioneering work in the 



early 80's treated the thick disk either as an outer flat- 
tened but exponential halo (fo r NGC 4452 and NGC 4762: 
lHamabe fc Wakamatsul Il989fl . or as a spheroidal bulge 
comp onent (for NG C 1381: Ide Carvalho fc da Costa , 
I l987h ffor NGC 3115: ICapaccioli. Held, fc NietoL 11987 : 
ISilva et a.lllT~ 



The detections of possible halo or thick disk stel- 
lar light in disk galaxies of later Hubble type have 
been made in a scant handful of mostly nearby 
edge-on galaxies (e.g. for ESO 342-017, IC5249, 
NG C 891, NGC 4565 NGC 5 907, and NGC 6504: 
der Kruit fc iwR ll 981 albt Tshaw fc GilmoreL FflRC . 



van 



Morrison. Boroson. 



- & Hardine 

1994; Sack ett et all 



)94t Iv an Dok knm et al 



Morrison et al 



Naslund fc JorsateiT Il997t iLeaueux et all 
Zheng et all Il999t iNeeser et all 120021) 



1997 



1998; 



Quite recently, 

Dalcanton fc Bernstein! l|2002t) suggest the detection of 
extended, ubiquitous thick disks in a large sample of late- 
type, edge-on galaxies by means of multi-colour imaging. 
However, their thick disks are solely detected by vertical 
colour gradients for which dust extinction complicates the 
interpretation. In addition, their vertical colour profiles, 
especially (R — K) , typically extend out to only very few 
vertical disk scaleheights. At those z-heights where they 
attribute the red colour to an additional component the 
thin disk may still be dominant over a potential thick 
disk and even determine the measured colour. 

In this paper, we analyse a set of eight edge-on SO 
galaxies using the classical approach for the identification 
of thick disks in external galaxies: the need for an addi- 
tional disk component when attempting to fit single disk 
models to the light distribution in a deep image. Thereby 
we characterise the presence and properties (scalelengths, 
-heights, and central surface brightnesses) of thick disks 
by directly measuring their structure. 
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2. Sample Selection & Observations 

The data were taken as part of the PhD s tudy on the ra - 
dial structure on galactic stellar disks bv iPohlenl (2001). 
The galaxies were selected according to the allocated ob- 
serving time, observatory, and CCD-chip size meeting the 
following morphological selection criteria. Using images 
from the Digitized Sky Survey (DSS) we verified that 
they are edge-on (i 86° ) , undisturbed, and similar to 
some disk-prototypical cases like NGC 4565 or IC 2531 to 
make it possible to consistently fit the applied simple disk 
model (cf. Sect.|3J). Galaxies with the following characteris- 
tics were rejected: spiral arms (indicating towards a lower 
inclination), a significantly asymmetric or disturbed disk 
(indicating towards strong interaction), and two-sided or 
significantly one-sided warped disks. In addition, galaxies 
which seemed to be dominated by the light of the bulge 
component, or appeared to be too small (-D25 ~2"), or 
showed only a faint, patchy disk were also excluded. For 
the observed lenticulars we therefore selected galaxies with 
a symmetric, smooth disk and a distinct, non-dominant 
bulge component. It is not unexpected that half of these 
galaxies show a box- or peanut-shap ed (b/p) bulge com- 
ponent, since iLiitticke et alJ l)2000a|) find that > 40% of 
all bulges are b/p shaped. Global properties of the finally 
observed eight SO galaxies are given in Tabled We do not 
argue that this sample is fully representative of the general 
population of SO galaxies but it ensures the best prospects 
for obtaining consistent models with our 3D modelling 
technique for all galaxies. 

The images (in Johnson R or V filter) were obtained 
in four observing runs in 1998/1999, three at the Danish 
1.54m telescope of the European Southern Observatory 
(ESO, Chile) and one at the 1.23m telescope on Calar 
Alto (CAHA, Spain). During all three runs at the ESO the 
1.54m Danish telescope was equipped with DFOSC and 
the C1W7/CCD which is a 2k x 2k LORAL chip providing 
a field size of w 13' and a scale of « 0.39" pixel -1 . The 
run at the Calar Alto 1.23 m telescope was done in service 
mode with the Site#18b chip, a 2048x2048 SITE CCD 
with 24 /im pixel size, providing an unvignetted field of 
?» 10' and a scale of « 0.5" pixel -1 . 

The standard CCD reduction techniques — overscan 
correction; subtraction of remaining large scale gradient 
in combined, oversan-subtracted, masterbias image; and 
careful flatfielding- were applied using the IRAF data re- 
duction package. Neither the DFOSC nor the Calar Alto 
CCD R-band images were affected by fringing. The in- 
dividual, dithered, reduced short exposures (150s-600s) 
were combined to the final deep image using IRAF's im- 
combinc task. These images are rotated to the major axis 
using the smallest angle of rotation according to their true 
position on the sky. Table|21summarises the detailed obser- 
vational para meters. During t he two ESO observing runs 
in 1999 several lLandold l)l992|) fields, partly enriched with 
additional stars provided by B. Skiff (priv. comm.), were 
observed. The standard fields were taken at least three 
times a night at different airmasses to determine the at- 



mospheric extinction. During the other two observing runs 
no standard stars were taken. The ESO run in 1998 is cal- 
ibrated by literature values and for galaxies without cata- 
logued values interpolated according to the measured sky 
background. Only a rough zero point could be estimated 
for the Calar Alto run by comparing a galaxy also ob- 
served in another calibrated Calar Alto run. For more de- 
tails a bout the photometric calibration we refer to lPohlenl 

3. Extraction of the Disk Parameters 

3.1. 3-Dimensional Disk Model 

We have developed a semi-automatic recipe to fit true 3- 
D single-component luminosity distributions to the 2-D 
data of edge-on galaxies and determine the galaxy pa- 
rameters, such as scalelength and scaleheight in a physi- 
cally meaningful way . Ou r method is de scribed in detail in 
IPohlen et alTl(2000bh and lPohlenl l|200lh and is only briefly 
recalled he re. The disk model is based upon the fundamen- 
tal work of Ivan der Kruit fc Searlel 1 1981aL They tried to 
find a fitting function for the three-dimensional light dis- 
tribution in disks of edge-on galaxies using the empirically 
determined exponential radial gradient, /ocexp(i?), and 
adding a description for the vertical distribution, /(z), of 
the stars. The luminosity density distribution L(R, z) can 
be written as: 



L(R 7 z) = L exp ( -— ) f n (z, Zq) H(i? co - R) 



(1) 



L being the luminosity density in units of [Lq pc ], 
Lq the central luminosity density, R and z are the radial 
resp. vertical axes in cylinder coordinates, h is the radial 
scalelength and zq the scaleheight, and n the index of the 
vertical distribution function. H(a;o — x) is the Heaviside 
function and R co is the cut-off radi us characterising the 
observed outer rad ial truncations ijvan der KruitL Il979t 
IPohlen et all I2OO2I and references therein). To limit the 
choice of parameters we restrict our models to the three 
main density la ws for the z-di s tribut ion (exp, sech, and 
sech 2 ) following Ivan der Kruitl l|l988lL Due to the choice 
of our normalised isothermal case zo is equal to 2h z , where 
h z is the usual exponential vertical scale height preferred 
by many authors: 



h{z) 

Hz) 



4 exp -2 - — '- 

zq 

2 sech I — 

,z 



sech ( — 
zo 



Fo r any details a bout th e numerical realisation we refer 
to IPohlen et all l|2000b|) . Therefore six free parameters 
(i, n, Lq, R co , h, zo) At the observed surface intensity on 
the CCD chip to the model. This model assumes that the 
vertical distribution is independent of position along the 
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Table 2. Observing log for the individual combined images with (i) the galaxy, (&) the filter, (s) the observing date, 
(4) the site, (5) the total coadded on-source exposure time i; nt , (5) mean seeing conditions during the observations, (7) 
the number of individual images with their ti nt . 



major a xis as is known to be true in good approximation 
(cf. e .g. Ivan der Kruit fc SearleL Il981at IShaw fc Gilmorel 
Il990|) . The increase of scaleheight with gal a ctoce ntric dis- 
tance as reported by Ide Griis fc Peletierl (|1997l) can be 
described as a combination of two disks each with con- 
stant scaleheight but different scalelength (cf. Sect. 13. 2} . 
Any deviation from constant scaleheights should be visible 
in the vertical profiles overplotted by the models, which is 
not the case (cf. Appendix[5J|. The apparent small vertical 
shift of some models in the lowest plotted vertical profile 
is due to the change in radial scalelength, as described in 
Sect. HO 

The possible influence on the six free parameters of the 
dust distribu tion, which was negl ected during the fit, is 
estimated in lPohlen et all l|2000b|) . There we have shown 
that even for a worst case scenario (large optical depth 
tr and a radially and vertically fairly extended dust lane) 
our model is able to reproduce the input parameters with 
an error of a typical 20%. We expect this effect to be 
even less significant for the present sample of lenticular, 
dust-depleted galaxies. Deriving individual errors on all 
parameters is a complex task since the main source of 
error is not the numerical fitting procedure (< 1%) but 
the systematic uncertainties in the process of fitting a 
rather simple, empirical model to real galaxies . An es- 
timation of the errors is given in IPohlen et all l)2000b(l 
by changing the applied boundaries of the region used to 
fit the data to the model (cf. also IPohlen. 2001). They 
found differences in h and zq of about 15% and Lq varied 
ab out a factor of two. This is in the same range found 
bv lKnaoen fc van der Kruiti l|199l|) when they compared 
published values of scalelength measurements from differ- 
ent studies. 



3.2. Thin /thick disks 

IPohlen! ll200lh noted that the SO galaxies in his sample are 
not well described by any combination of a single disk and 
another spheroidal component, such as a de Vaucouleurs 
R 1/A bulge model l|de VaucouleursL Il94sh . All SO galaxies 



reveal a typical, continuous change of slope when one com- 
pares the major axis with parallel profiles above/below the 
plane. The profiles significantly flatten towards cuts higher 
above the midplane. According to the single component 
model L(R, z) all slopes should be nearly parallel. We in- 
fer from the images (cf. Appendix [5J) that all SO galax- 
ies show a kind of smooth outer envelope or highly flat- 
tened spheroidal component. This deviation from a nor- 
mal shape cannot be explained by a bulge component. Any 
large i? 1 / 4 bulge would have to be apparent on the major 
axis which is certainly dominated by an exponential (disk) 
component. It is worth mentioning that in the sample of 
IPohlenl l|200ll) there are four galaxies, NGC 3390 (T = 3), 
NGC 3717 (T = 3), NGC4696C (T = 3.4), and NGC 6504 
(T = 2) a that showed a similar behaviour but were classi- 
fied as late-type galaxies. 

The outer com ponent could be described as a 'thick 
disk' according to iBursteinl l)l979cj) with a flatter slope 
than the inner thin disk, equivalent to a larger scalelength. 
However, the main characteristic of a thick disk is that the 
observed vertical profiles depart from the simple exponen- 
tial model. Compared to the model, the outer parts are 
systematically brighter wit h increasing distanc e to the ma- 
jor axis. As mentioned bv lde Griis fc Peletierl (^997), the 
measured increase of the scaleheight with radius ('flaring') 
for early-type, edge-on disk galaxies can be understood if 
these galaxies have both thick disk scalelengths and scale- 
heights larger than for the dominant old disk. 



3.3. Fitting Method 

To apply our single component fitting method described 
in the previous section we have to assume that there is a 
region in the galaxy dominated by only one of the com- 
ponents. If one uses exponential profiles it is obvious from 
the vertical and radial cuts that there are well separated 



a The classification i s done by |Lutt4cke_ei L _aI . (2000a). It is 
consistent with that of Ivan Dokkum et^hTEfflffi. but signif i- 
cantly different from that of the UGC catalogue TNilsonllT973Tl . 
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vertical ranges (around the major axis and the outermost 
profiles) in which the light of one of the two disks dom- 
inates. Fitting thin and thick disk component simultane- 
ously would be the desired approach. However, our sin- 
gle comp onent model ha s already six free parameters. As 
shown in Pohlen l)200l[) . fitting this model to observed 
data is a non-trivial task. The problem is the applica- 
tion of an idealized model itself, which obviously will not 
totally accurately describe the measured two-dimensional 
light distribution. This technique requires continuous hu- 
man supervision to control the influence and quality of 
each individual parameter. This would most probably be 
even more difficult for a parameter set twice as large. We 
decided therefore to use an iterative fitting routine start- 
ing from outside-in since the thick disk clearly dominates 
the outer parts. The first step is to determine an initial 
estimate for the outer disk by restricting the region to be 
fitted and using our single component model. The next 
step is subtracting the derived full thick disk model from 
the original image and fit an inner disk to the residual 
by restricting the fitting range to the inner parts. Then 
we start from the beginning by subtracting this inner disk 
model from the original image and fit again the outer disk, 
now to this residual image. The initially pre-defined fitting 
regions for the thin and thick disk are sometimes adapted 
slightly after subtracting one of the components and be- 
fore starting the second fitting round. As it turned out 
this process is remarkably stable for most of the galaxies. 
After one iteration the disk parameters are already the 
same within the range of one single model. The reason 
for this is the domination of the thick disk at large ver- 
tical z-heights ( J5 3.5zq = measured here with high 
S/N, and the large radial range used to fit the thick disk 
component. 

To restrict the number of free parameters we decided 
to prefix the inclination during the fitting process for this 
sample. For most of the galaxies the symmetric shape of 
the bulge and disk component indicates an exactly edge- 
on orientation. Only one of the SO galaxies (NGC 3957) 
exhibits a dust-lane making it possible to estimate the 
inclination at i = 88° ± 1° following the method by 
iRarteldrees fc Dettmarl fiflfl^. 

The different density laws for the vertical dis- 
tribution are similar for large z and only differ 
around the mid-plane of the luminosity distribu tion 
(cf. Fig.4 inlde Griis Peletier. fc van der Kruitl .lT99^. In 
Pohlc n et al. ( 2000b) the actual choice of one fitting func- 
tion is done individually for each galaxy depending on 
the measured profiles. However, near the plane of the 
galaxy the contribution of a thick disk is much smaller 
than that of to the thin disk component. First tests with 
a free choice of the fitting function for the vertical den- 
sity distribution (exp, sech, or sech 2 ) showed in the case 
of ESO 311-012 that the iterative fit became unstable for 
this reason. In the subsequent modelling we have chosen 
the intermediate fitting function f2{z) ocsech(2z/zo) in all 
cases for the thin and thick disk component. There was 
only one galaxy, NGC 2310, for which the sech function 



did not yield a satisfying convergency of the iterative fit- 
ting. However, switching to sech 2 significantly improved 
the fit. Finally, only three free parameters are left for each 
disk (L , h, z ). 

We want to point out that these decompositions are 
thus model dependent. An isothermal thick disk con- 
tributes less light to the thin/thick disk combination than 
to an exponential one and the shape of the thin disk will 
also be different. In addition, vertical scaleheights zq ob- 
tained with an exp model are systematically larger than 
in a sech model and again larger than in a sech 2 fit. 
Depending on the vertical boundaries chosen for the fitting 
these differences could be more than 30%. The choice for 
the vertical distribution influences also t h e bes t-fit scale 
parameter ratio z^/zq ■ IShaw 

fc Oilmorel ifmsi used all 
combinations of either sech 2 or exp models for NGC 4565 
and derived Zq/zq in the range of 4.3 — 5.4. 

In contrast with some of the previous 1-D-only fitting 
methods we do not simply use the deviation from a simple 
exponential on the minor axis, or parallel profiles; to mea- 
sure thick disks we use the full 2-dimensional information 
and are therefore able to fit the radial scalelength of the 
thick disks. 

3.4. Breaks in Radial Profiles 

One of the main difficulties while properly fitting our 
model to the SO galaxies is their outer disk structure . 
This will be discussed in detail in IPohlen et al.l l|2004h 
and only briefly addressed here. As one can see from 
the figures in Appendix [5] there are clear breaks in 
the outer parts of the radial profiles. These are simi- 
lar to the tr u ncatio n s in more later-type galaxies (cf . 
[Pohlen et all Eoo3 iKretrel. va,n der Kruit. fc de GriisL 
h(M Ide Griis. Kregel. fc Wessorl Hflml 

IPohlen. Dettmar. fc LiittickeL 12000^ - For all galax- 
ies (except ESO 311-012) a similar break, slightly less 
pronounced, is also seen in other radial cuts parallel to 
the major axis. Our model, according to Eq. Q), however, 
describes only an infini tely sharp truncation R co . As 
shown in IPohlenl l)200lh this implies a tight coupling 
between the radial scalelength h and the cut-off radius 
-R co , when using our model fitting data with the observed 
breaks in the profiles. In addition, the sharply truncated 
model exhibits an i ntrinsic bendin g of the profile towards 
the outer parts (cl. IPohlenl l200lh . This complicates the 
visual quality control compared to the more flat infinite 
exponential model without any truncation. Therefore we 
decided to use the infinite exponential model, realised 
within the same fitting program by fixing R co to ten 
times the radial scalelength. For our thick/thin fitting we 
restricted the fitting region to points inside the observed 
break radius. However, fitting an infinite exponential 
model to the intrinsic two- slope p r ofiles is also affected 
by s ystematic errors (cf. IPohlenl l200lt IPohlen et all 
Depending on the ratio of the inner, shallow slope 
(h- m ) up to the break, to the steeper, outer slope (/i ou t ) 
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beyond the break radius the determined scalclcngth will 
be systematically too small. Assuming a mean ratio of 
/tin/ftmi t. = 4.4 ±1.7 from a large, edge-on sample by 
Pohlen (200l| makes it possible to quantify the expected 
offset. The best fitting scalelength /i; n will be about 26% 
larger compared to the intrinsic radial scalelength (cf. 
IPohlen et all 12004) . 

Although the exact value of the scalelength is thus 
model dependent, neither zq nor /in are influenced by this 
problem. The derived scaleheights are independent of the 
scalelength and therefore robust, but of course still de- 
pending on the chosen density law (exp, sech, or sech 2 ) 
for the z-distribution. 



4. Results 

For six (75%) out of the eight SO galaxies we were able to 
derive consistent thick/thin disk solutions (cf. Sect. 14. 
The resulting best fit parameters are listed in Table 
and Table 0] for the thick and thin disk, respec- 
tively. The resulting radial and vertical profiles overplot- 
ted by our best fitting model are shown in Appendix 
IbI The derived thin disk scalelength and scaleheight 
valu es are consistent with those found in studies 
bv iKregel. van der Kruit. fc de Griisl ll2002l) and IPohlerJ 
l|200lf) . containing mostly galaxies of Hubble type later 
than Sa. 

The mean ratio of thick to thin scaleheight for the 
five galaxies with the same vertical model combinations 
(sech+sech) for the thick and thin disk is Zq/zq = 3.6±1.0 
and for the scalelength h k /h n — 1.8 ± 0.1. Including also 
NGC 2310 we find Zq / Zq = 3.4±1.0 and h k /h n = 1.9±0.4. 
We derive central surface brightnesses (uncorrected for in- 
clination) in the range of 22.0 < fx k < 22.2 R-mag/n" and 
21.9^4 £ 22.9 V-mag/n" for our thick disks compared 
to 19.2 </xg < 19.6 R-mag/n" and 19.5 <^ ^20.7V- 
mag/n"for the thin disks. This implies that the contri- 
bution of the thick disk to the central surface brightness 
is about 10% of that of the thin disk. The thick disk cen- 
tral luminosity density Lq ranges between 3.5% and about 
10% (mean: 5.6%) of the thin disk value. The ratio of the 
total luminosities of the thick and thin disk is between 
about one third and one. 

The profile on the minor axis in addition to radial cuts 
high above the major axis (cf. Appendix[5J) reveals that for 
four (NGC 2310, NGC 3564, NGC 4179, NGC 5047) out of 
the six fitted galaxies there is no significant bulge compo- 
nent visible at large vertical height above the disk. Even 
the apparently bulge dominated galaxy NGC 3564 is well 
described by a thick disk component. The remaining bulge 
light after subtraction of the thin/thick disk combination 
is comparable to that of the other galaxies. This implies 
that all bulges of our SO galaxies could not be well de- 
scribed by a traditional de Vaucouleurs i? 1 / 4 bulge. In ad- 
dition, any bulge component would be too flat to account 
for all the light high above the disk at large galactocentric 
radii. 
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Table 5. Radial and vertical fitting regions for the 
thick/thin disk components with (1) galaxy name, (2) be- 
ginning and end of thick/thin disk vertical region, and (s) 
beginning and end of thick/thin disk radial region 



4.1. Fitting regions 

One of the important constraints for fitting empirical, sur- 
face photometric models to observed data is the definition 
of the actual regions which are marked to characterise the 
individual model components. Therefore we list in Table 
0for each galaxy the radial and vertical ranges where our 
thick/thin disk fit was applied. In all cases we restricted 
the two fitting regions to be distinct from each other and 
outside the dust lane (in the case of NGC 3957) or any 
bar/ring like feature visible in the radial profile. This is 
obvious for the first thick disk estimation. However, in 
principle one could extend the fitting range towards the 
inner/outer parts for the following thick/thin disk itera- 
tions, since the other component in each case has already 
been subtracted. In addition, we masked by hand bright 
stars and background galaxies, forcing the program to ig- 
nore these regions. 

4.2. Comparison with Literature 

iBursteir] l|l979oj) only describes the 'thick disk' component 
qualitatively as being more diffuse than that of to the 
inner, dominating disk, and possessing a flattened shape. 
The important observation is that the scaleheight of a 
fitted exponential to the vertical profiles increases with 
radial distance R. He determines a flattening, defined as 
the intrinsic z-thickness compared to the diameter a, of 
z/a = 0.25 — 0.2 at /xb=25.0 which is somewhere between 
thin disks (z/a 0.1) and E4 ellipticals (z/a = 0.6) for five 
edge-on SO galaxies. We find values of 0.29 (NGC 2310) or 
0.32 (NGC 5047) comparing minor to major axis diameter 
at the second outermost contour (cf. Fig. lB.2IB.8|) . 

While comparing scale parameter ratios in the liter- 
ature one has to keep in mind that they are often ob- 
tained with very different fitting methods and even dif- 
ferent fitting functions for the vertical distribution of 
thin and thick disk light. As discussed in Sect. 13.31 this 
already implies systematic differences of at least 20%. 
Surprisingly, there are no detailed parameter studies for 
thick disks of SO galaxies in the literature that provide 
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Table 3. Results for the thick disk: (i) Galaxy, (2) distance, (3) filter, (4) central surface brightness (uncorrected for 
inclination) of the thick disk model [1q , (5) vertical scaleheight Zq in arcsec, parsec, and in units of the thin disk vertical 
scaleheight Zq , and the (6) radial scalelength /i k in arcsec, kpc, in units of the thin disk radial scalelength h n , and in 
units of the thick disk vertical scaleheight Zq a (used sech 2 instead of sech) 
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Table 4. Results for the thin disk: (1) Galaxy, (2) distance, (3) filter, (4) central surface brightness (uncorrected for 
inclination) of the thin disk model fj,Q , (5) vertical scaleheight Zq in arcsec, and parsec, and the (6) radial scalelength 
/i n in arcsec, kpc, and in units of the thin disk vertical scaleheight Zq 



both scaleheight and scalelength ratios (cf. Sect.QJ. The 
only scaleheight ratios for SO galaxies (cf. Table |HJ) are 
giv en for two galaxies (NGC 4710, NGC 4762) in detail 
bv Ide Griis fc van der KruiH (^996) and for five galaxies 
with Hubble type TJSO.O without prov iding the individual 
values or applied fitting functions bv lde Griis fc Peletierl 
l|l997(l . They find ratios of z$/ 'z% = 1.8 - 4.6, which falls 
well within the range of values for our SO galaxies. 

There are clearly more studies in the literature that 
provide detailed parameters (scaleheights and often also 
scalelengths) for thick disks in galaxies of later morpho- 
log ical type. 

IVan Dokkum et ahl l|l994|) modeled excess light at 
large vertical distance one-dimensionally with a thin/thick 
disk combination for the edge-on Sab galaxy NGC 6504. 
They found that the ratio of scaleheights is roughly 4, 
slightly higher than the mean value for our SO sample. 
However, their central surface brightness for the thick disk 
is significantly lower (A/zo = 1.6 mag) than our value for 
NGC 3957. This effect is not related to the different fit- 
ting functions used since our sech model should tend to 
result in systematically lower values for /iq than in their 
exp model. 

For the Milky Way, recent optical star-count measure- 



ments bv lLarsen fc Humphreys! l)2003|) yield thin disk val- 
ues of Zq — 300 pc and h n = 3.5 kpc, and for the thick 
disk Zq = 900 pc and h k = 4.7 kpc. These are s imilar 
to the infrared 2MASS star counts by lOihal l|2001^ with 
zj =260pc, h n =2.8 kpc, z% =860pc, and h k =3.7 kpc, 
keeping in mind that the near-infrared surveys always de- 
rive systematically smaller scalelengths than in the opti- 
cal. The resulting ratio of thick to thin disk scaleheight is 
very similar to the mean value we find for our S O galaxies 
but th eir ratio of the scalelengths is also smaller. IPu et all 
l|2003|) present a list of thick disk local normalizations (rel- 
ative to the solar neighbourhood) obtained in different 
studies including their own measurements (cf. their Table 
2). The quoted values range from 2% up to 13% with a 
mean value of about 6.1%, which agrees very well with the 
mean value of 5.6% for the central luminosity density of 
ou r SO galaxies. 

iMorrison et al ] lll997h used one-dimensional vertical 
fits to derive the thick and thin disk scaleheights of 
NGC 891, an edge-on Sb galaxy similar to our Milky 
Way. They quote a large range of possible Zq/zq values 
which are consistent with our SO d i sks. P revious mod- 
eling of NGC 891 by Ivan der Krvul (|l984 yield a ratio 
of Zq/zq = 3.0 for a three component (thin plus thick 
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Table 6. Comparison with literature: (i) Galaxy name (s) morphological type, (s) ratio of thick to thin disk scale- 
height, (4) ratio of thick to thin disk scalelength, (5) central surface brightness of thick disk (if quoted), (6) centra l 
surfa ce brightness of thin disk (i) refer ence DE96lde Griis fc van der Kruitl lll996t. DE 97 lde Griis fc Peletier] (1997), 
VA94 lvan Dokkum et all lll994h . WU f)2lWu et alJ ll2002t. MQ97lMorrison et alJ dl997h . VA84 lyan der Kruitl 1 1984 1. 
SH89 lShaw fc Gilmorel (Il989l) . LAOs lLarsen fc Humphreys! ll2003h . OJOl lOihal d200lh . NE02 lNeeser et all (l2002t . 
AB99 lAbe et all l|l999|) 

a NGC 2310 excluded (only sech 2 model for thick disk) c three component model (thin disk, thick disk, and halo) 
b three component model (bulge and thin plus thick disk) d starcount analysis 



disk, and additional i? 1 / 4 -spheroid) fit which is also in 
the range of our SO galaxies. The slightly higher value 
of the thick disk central luminosity density (17% of the 
thin disk) compared to our maximum of 10% could be 
explained by the additional contribution of the applied 
i? 1 / 4 -sphe roid, especially hig h abov e/below the midplane. 
However, iBah call & Kvlafis! 1^985) already pointed out 
that in the case of NGC 891 a combination of thin plus 
thick disk and a model with only one disk and a i? 1 / 4 - 
spheroid both give a good description of the data used by 
Ivan der Kruil l|l984|) . 

For another prominent edge-on Sb galaxy, NGC 4565, 
several attempts for a de composition are available. More 
recently. I Wu et all 1 20021) presented a thin/thick plus halo 
decomposition from a deep intermediate-band (6660A) 
image deriving values only slightly smaller than for our 
SO galaxies (cf. Table ■ However, the two comp onent 
disk models of NGC 4565 bv lShaw fc Gilmorel ill989f) yield 
significantly hig her values for the ratio Zq/zq . The light 
IWu et all l|2002j) put down to a halo is here ascribed solely 
to a thick disk with a large scaleheight. There is a large va- 
riety of mod el combinations (with different components) 
possi ble (cf. IShaw fc Gilmorel Il989t iNaslund fc Jorsateil 
Il997|) . which complicates a direct comparison. 

NGC 5907, another nearby edge-on Sc galaxy, is 
an even more complex case. It is not yet clear 
if the multiple detected extended light distribution 



is a thick disk or a stellar ha l o component (cf . 
Mor rison. Boroson. fc Hardinglll99l . lZheng et alJ (Il99fll) 

even concluded that NGC 5907 does not have a faint 
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obtained a deep optical image of the 
edge-on Sd galaxy IC 5249. They detected additional light 
to that predicted by a single exponential disk and tried to 
fit a thick disk model. The scaleheight ratio of their best 
thick disk is well within the range for our SO galaxies. 
However, the scalelength of the thick disk is exceptionally 
smaller compared to the thin disk and their central lumi- 
nosity density (Lq/Lq = 7.4) is noticeable larger than for 
our sample. 

Recently, iNeeser et al.1 l)2002l) reported for the first 
time the detection of a thick disk in a low surface bright- 
ness galaxy (ESO 342-017 classified as Scd). They find for 
the scaleheight a ratio of Zq/zq =2.5 (close to our val- 
ues for NGC 2310 and ESO 311-012) with a comparable 
or somewhat larger scalelength for the thick disk. The 
scaleheight ratio and unprojected central surface bright- 
ness (cf. Table [SJ are surprisingly similar to the range of 
values found for our SO galaxies. 
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5. Discussion 

5.1. Thick disks: discrete or continuous? 

A key question is: How to describe a thick disk in gen- 
eral? With reference to the proposed different formation 
scenarios described in Sect. ^ we can assume that the 
disk component is either characterised as a superposition 
of two discrete and independent isothermal disk systems 
(as done here), or built from the con tribution of multiple 
velocity-dispersion comp onents (e.g. IWielen et all Il992l 
iDove fc Thronsorlll993(l . These two approaches on how to 
treat a thick disk seem to be incompatible. Therefore any 
derived parameters for the two-disks model appear useless 
for the multi-component disk. In addition, the latter seems 
to be superior since it is consistent with the model of con- 
tinuous disk heating leading to t he well known observable 
age- velocity dispersion relation JWielenlll977h . 

However, M aiewskil (|l993j) already states in his re- 
view that detailed studies of the spatial and kinematical 
distribution of stars in the Milky Way do not make it 
possible to decide if the thick disk is a discrete compo- 
nent or a more continuo us sequence of stellar populations. 
Rece ntly. iNissenl ll2003l) concluded that the latest studies 
(e.g. iBensbv. Feltzing. fc Lundstroml l2003h argue for the 
separate entity picture. This confirms our result that 75% 
of the chosen SO galaxies are well described with a distinct, 
two-component thin/thick disk system. Note that these 
results exclude all thick disk formation scenarios based 
solely on heating . Especially the elemental abundance 
trends found by IBensbv. Feltzing. fc Lundstroml ((2003) 
favour a merger scenario where a satellite galaxy ei- 
ther merges with the parent galaxy or sheds significant 
a mounts of its materia l to form the thick disk as proposed 
bv lAbadi et aD l|2003h . 

However, as already noted by Ijacobi fc Keeell (|l994^ 
it does not seem possible to distinguish between our sim- 
ple model and more sophisticated ones. Therefore the very 
good description with the applied discrete two-component 
system for our galaxies is only a first step towards study- 
ing the luminosity distribution of external galaxies. In ad- 
dition, one has to keep in mind that our model only fits 
well in a restricted radial range (only out to 2.2h n — 3.9h n ) 
where most of the inner part is "hidden" by the central 
bulge, bar, or ring components. One key to the nature of 
the thick disk may lie in studying the outer parts where 
the breaks in the radial profiles of both thin and thick 
disk should provide an additional constraint. Any vertical 
colour gradients could provide additional information but 
the only avail able colour map for our sample (NGC 3957 in 
IPohlerl EoOljl suffers from low S/N and is not conclusive. 

5.2. Are thick disks of SOs exceptional ? 

The fact that the range of thick disk parameters for all 
known SO galaxies are not too different from those of late- 
type galaxies, even compared to a low surface brightness 
Scd galaxy, is especially surprising, since at first glance 



one expects SO galaxies to posses more prominent thick 
disks. However, in terms of the thick-to-thin-disk scale- 
height ratio our values agree well with those derived for 
all other galaxies (cf. Table EJ). Even taking into account 
the central surface brightness or central luminosity den- 
sity the comparison yields similar values. Does this point 
to a general formation process for thick disks independent 
of bulge-to-disk ratio and Hubble type ? At this stage we 
are not able to answer this question. Although the num- 
bers of galaxies used for the comparison are the same, all 
literature values come from different sources, using some- 
times very different methods and models to derive their 
parameters. In the case of NGC 4565 it is obvious that de- 
pending on the model, one can find an even larger variety 
of Zq/zq -values than the full r ange for our six galaxies. 
In ad dition, as pointed out by iKnaoen fc van der Kruitl 
and again discussed in detail bv iPohlenl (|2001), 
fitting disk-model components to surface brightness data 
of edge-on galaxies is a delicate business, and comparing 
the results of different authors can be misleading. To over- 
come this problem one has to extend this survey of thick 
disk parameters also to late-type galaxies. To reduce the 
inevitable influence of their prominent dust lanes this has 
to be done in the near-infrared. After applying the same 
fitting method one is able to address the questions if thick 
disks are really common around all Hubble-types, and if 
their parameters are really similar, as suggested here. This 
would entail a common formation scenario for thick disks 
independent of their normal evolution along the Hubble 
sequence. 

Of our galaxies with successful thick disk models there 
are only two for which there are rotational velocity mea- 
surements in the literature. Therefore any correlation be- 
tween mass and thick disk parameters could not be de- 
rived. Consequently it is still unclear if the mass, as a 
galaxy characteristic, is related to the thick disk parame- 
ters. 

5.3. Thick disk scalelengths 

We find significantly larger scalelength ratios (h k /h n = 
1.9) for the thick disks of our SO galaxies than in the 
literature (h k /h n ^1.4). Is this a possible distinction in 
thick disk parameters between late-type and our early- 
type galaxies ? Again, differences in fitting methods, espe- 
cially for the scalelength in edge-on galaxies, could be re- 
sponsible for this apparent disagreement. In particular, fit- 
ting profiles with clear breaks using a single disk (with an 
infinitely sharp truncation or infinite exponential as done 
here) will entail systematic errors (cf. Sect. 13.4(1 . therefore 
this question could only be addressed unambiguously ap- 
plying exactly the same method to late-type galaxies. 

However, although the exact number is uncertain the 
scalelength of the thick disk is without doubt systemat- 
ically larger than that of the thin disk. Except in one 
case this is also true for all the literature values. Does 
this larger scalelength imply a different formation pro- 
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cess or even restrict the available formation scenarios ? 
At a first glance the different scalelengths contradict a dy- 
namical heating scenario, since this should only alter the 
vertical distribution. However, assuming two distinct disk 
components with radially cons tant thickness the velocity 
dispe rsion scales with radius (|van der Kruit fe Freemanl 
Il986lh It depends therefore on the exact way (radial dis- 
tribution) the proposed mechanisms (cf. Sect.Pl dispenses 
the energy in the vertical motions of the stars. In the 
case of heatin g the disk by satellite accretion , the N -body 
simulations of iQuinn. Hernauist. fc Fullaearl l|1993|) show 
that the scalelength of the disk is nearly unchanged (even 
slightly smaller than the initial scalelength) for the inner 
parts (out to 3/i ini ). Due to the migration of material out- 
wards in radius, at larger radius (out to ~6/ij n i) the final 
disk sho ws a second shal l ower c omponent (cf. their Fig. 4). 
However. lAeuerri et alJ <)200l[) find in their N-body simu- 
lation a global outward transport of disk material leading 
to a general increase of the disk scalelength of 10% to 60%. 
Although this increase is below our measured scalelength 
ratios, the latter seems to coincide in general with observa- 
tions of larger thick disk scalelengths. One has to keep in 
mind that in the satellite accretion scenario a pre-existing 
thin disk gets heated. The thin disk to be observed to- 
day must be rebuilt out of gas remaining after the merger 
process. It is not guaranteed that this thin disk will have 
the same scalelength as its predecessor since we do not 
have a unique explanation of how the disk scalelength is 
determined out of an initial (or new) gas distribution (cf. 
IPohleu et a,lll2000hh . 

Looking from a slightly different angle, one can also 
try to apply the larger thick disk scalelength as a general 
argument against an internal heating scenario. Assuming 
infinitely exponential disks, the larger scalelength implies 
a larger angular momentum content of the thick disk. Any 
valid heating mechanism must therefore add angular mo- 
mentum, which is not the case for the internal heating 
scenarios. However, note that the disks are not at all in- 
finitely exponential but show clear breaks in their profiles 
(cf. Sect. IH.4jl . The disks exhibit the break at roughly simi- 
lar radial radius (e.g. Fig. lB.8|l . In this sense the thick disk 
is truncated "earlier" (in respect to its scalelength) than 
the thin disk. In addition to a probably lower rotational 
velocity this could add up again to a similar angular mo- 
mentum content for thick and thin disk, ruling out the 
whole argument against the heating scenarios. The key 
point here is again the origin of these breaks in the radial 
profiles. 

5.4. Residuals 

In addition to comparing thick disks across a large range 
of characteristic parameters of galactic disks, fitting and 
subtracting a combined disk yields valuable information 
on the structure and size of galactic bulges. As discussed 
in Appendix ^ the residual images highlight the devia- 
tion of the galaxy from the disk model. Structures only 



faintly indicated in the profiles become obvious. One im- 
portant result related to this is that for our sample of 
SO galaxies none of the present bulges could be described 
as a traditional earl y-type R 1 ' 4 bulge in agreement with 
iBalcells et al J l|2003h . However, a detailed structural anal- 
ysis of the remaining central structures has to be done 
with great care. As discussed in Sect. 13.31 choosing the 
vertical model is not a unique process and therefore any 
variation of f(z) for the thick and thin disk could alter 
the shape of the disk profile in the inner bulge parts. 
Especially for galaxies where the bulge component along 
the minor axis extends vertically above the disk (ESO 311- 
012 and NGC3957) this could be even more uncon- 
strained. 
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Appendix A: Comments on Individual Galaxies 



NGC 1596: (Fig.|BJJ) A 'normal' SO enclosed in a huge, 
spheroidal, outer envelope possibly triggered by interac- 
tion with the nearby physical companion (NGC 1602). The 
difference from other SOs is obvious from radial cuts. It is 
the only galaxy in which the profiles bend upwards (not 
caused by a background gradient) at larger galactocen- 
tric distances (i? = ±115"). This looks similar to an R 1 / 4 
bulge component beginning to dominate the profile when 
extrapolating the curvature of the inner bulge towards the 
outer extension. However, the overall flatness of such an 
R 1 / 4 bulge would be unusual. In addition, the contours 
start to show substructure in the very outer parts (in con- 
trast to the smooth contours of ellipticals). The inner light 
profiles are dominated by bulge and a possible inner disk 
component. The major axis shows several distinct parts: 
an inner bulge, a steep straight decline out to i?«±80", 
an asymmetric slightly shallower part out to R= — 110" 
and R = +120", followed by an upwards bending profile 
an upbending profile. The SW-side is slightly more ex- 
tended than the NE-side. After the outer thick disk has 
been subtracted, the remaining profiles clearly resemble 
still the thin/thick disk behaviour (changing of slopes to- 
wards outer/higher radial cuts), therefore no consistent 
model could be derived. 

NGC 2310: fFig. lB~2ll exhibit s a box- shaped bulge (type 
2 according to iLiitticke et all l2000al) . The bulge/bar is 
rotated «0.5° against the outer disk, therefore we find 
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a slightly tilted residual. The radial profiles show a com- 
plex structure with inner, small, very bright nuclear bulge 
component, and further out a significant dip in the pro- 
file (caused by a possible dust ring or mimicked by a 
subsequent stellar ring), sitting on top of a bar/ring-like 
structure R ps ±50". There is a narrow disk- like struc- 
ture i?«±82", here assumed to be the main underlying 
disk, finally showing a straight decline into the noise (NE- 
side disturbed by nearby star-halo). The best fitting model 
overestimates the inner disk component due to inner ring- 
like feature. This feature is responsible for the apparent 
holes in the residual image, which shows the inner box- 
shaped bulge in addition to a large ring or bar. 

ESO 311- 012: fFigTOl Box-sh aped bulge (type 2 ac- 
cording to iLiitticke et all l2000aj) . The galaxy position is 
close to the Galactic plane with many foreground stars 
overimposed. The bulge component is still visible on the 
minor axis high above the plane. The radial profile shows a 
bar-like structure at i? = ±20", followed by an exponential 
disk component out to a break at i?«±83". The residual 
is a box-shaped bulge with the indication of a weak bar 
extending on both sides. 

NGC3564: CFig. lB~4)l The galaxy is probably part of 
small group/cluster with at least two large, nearby phys- 
ical companions (NGC3568 and NGC3557). Taking into 
account the low surface brightness part it is different from 
the other SOs. Out to /iy = 23.4mag/n"it looks similar 
although with a more dominating bulge. The outer con- 
tours, on the other hand, reveal that NGC 3564 is sitting 
in a huge, roughly spherical stellar envelope, which has 
probably been built up by a merger with a smaller satel- 
lite galaxy still visible in this structure on the NE-side 
((R, z) = -47", +24", distance unknown). The radial pro- 
file shows an inner bulge-dominated part, followed by a 
bar-like plateau i?«±10", sitting on top of an inner disk 
structure out to R « ±30"which is followed by a straight 
decline into noise. The residual resembles a rather large, 
spheroidal inner bulge together with a bar or ring-like 
structure extending to both sides. It nicely shows the pos- 
sibly merging satellite. The thick disk component takes 
most of the light from the envelope and therefore the re- 
maining bulge light is comparable to to that of the other 
galaxies. 

NGC 3957: fFig. lB~5)l Galaxy cluster in back- 
gr ound. Nearly box-sh aped bulge (type 3 according 
to lLiitticke et all l2000a|) . It is an SO galaxy with unusual, 
inner (R ± 30") dust lane. The bulge component is still 
visible on the minor axis high above the plane. The radial 
profiles show inner bulge, followed by a ring- or bar-like 
component at R w ±30" (visible in cuts above the dust 
lane) , with break at R ~ ±60" and followed by a straight 
decline into the noise. The profile probably starts to 
bend outwards from R ps ±105" (/iy ^24.3mag/n") to 
R sa ±150". The residual resembles a nearly box-shaped 
bulge with an indication of a bar extending to both sides. 
The apparent holes are caused by absorption in the inner 
dust lane. 



NGC 4179: fFig. lB~6|) The inner bar/bulge component 
is rotated «0.5° against the outer disk. The major axis 
shows a bright inner bulge, a disk- or bar-like feature be- 
tween RRi±40"(z<±5"), a more shoulder- like decline out 
to i?«±80", here taken to be the main underlying disk, 
followed by a straight decline into the noise. The resid- 
ual resembles a spheroidal bulge with an indication of an 
inner disk or bar extending slightly to both sides. 
NGC 4521: fFig. lB~7)l The inner component is ta 
1° rotated against the outer disk. The radial profile shows 
a bar- or ring-like structure between R « ±20"out to 
z ^ ± 5"; on top of a two-slope disk, inner part out to 
a break radius of 60", followed by a steeper straight 
decline out to R « — 77"and a slightly upwards bending 
part at the very edge (perhaps flatfield problem) . The only 
galaxy with inconsistent thick/thin disk fitting without 
obvious reason (cf. NGC 1596). All trials, such as changing 
vertical fitting function, have failed. When various thick 
disk models have been subtracted, the resulting thin disks 
always contributes too much light in the outer/higher z 
profiles dominated by the thick disk. Because of this the 
following thick disk is flattened, producing an unstable sit- 
uation, perhaps caused due to lower image quality, prox- 
imity to chip edge, and disturbed by bright star on the 
NW-side. In principle, the galaxy exhibits the same fea- 
tures in radial/ vertical cuts as our other SOs. 
NGC 5047: fFig. lB~8ll Pea n ut-shap ed bulge (type 1 
according to ILiitticke et all l2000aj) . There are many 
companions around associated with the foreground 
NGC 5044 group of galaxies (Av w 3600 kms" 1 ). Outer 
disk has a slightly s-shaped warp, on the SW-side 
towards N, on the NE-side towards S. The inner disk 
part is slightly off-centered to the outer component. The 
radial profiles show a bar-, or ring-like feature between 
i?« ±30"out to z ^ ± 3", on top of a shoulder-like inner 
disk out to R « ±50", followed by a straight decline 
out to fiy <; 25.0 mag/ □" . The residual shows a nice 
peanut-shaped bulge with a thin bar-like component. The 
ratio of the bar length to the maximum of the peanut 
distortion closely matches simulations explaining the 
peanut by dynamical processes in the bar potential (cf. 
ILiitticke et alll2000hh. 

Appendix B: Isophote maps, radial/vertical 
profiles, and residual images 

The following figures show in the upper panel an isophote 
map of each galaxy rotated to the major-axis. A small 
arrow in the upper right corner indicates the North. The 
magnitude of the outer contour (/Ltii m , defined by a 3<r cri- 
terion of the background) is indicated in each plot. The 
consecutive contours are equally spaced by 0.5 mag. The 
contour lines are plotted with increasing smoothing to- 
wards the outer parts. For the inner contours out to where 
the noise begins to increase no smoothing was applied. 
The following 2-3 contours are smoothed by replacing each 
pixel by the mean of 3 2 -pixels "around" the central pixel. 
For the outer two contours this smoothing is increased to 
5 2 -pixels. The second panel displays the major-axis sur- 
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face brightness profile (top) and two parallel radial pro- 
files each above and below the major axis. The exact 
z-positions for the plotted profiles are indicated in the 
upper left corner of the plot. For a consistent thin/thick 
disk decomposition the model profiles are overplotted as 
dashed lines. The third panel displays the minor-axis sur- 
face brightness profile (top) and two parallel vertical pro- 
files on either side of the minor-axis. Again the exact /im- 
positions for each cut are marked in the upper left corner 
and the best model is overplotted with dashed lines. The 
lower panel is the isophote map of the residual image af- 
ter subtracting the best fitting thin/thick disk model. The 
magnitude of the outer contour (fi\i m ) is again indicated 
in each plot and the following contours are equally spaced 
by 0.5 mag. 
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Fig. B.l. NGC 1596 R-band 
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Fig. B.2. NGC 2310 R-band 



Fig. B.3. ESQ 311-012 V-band 
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Fig. B.4. NGC 3564 V-band 



Fig. B.5. NGC 3957 V-band 
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Fig. B.6. NGC 4179 V-band 



Fig. B.7. NGC 4521 R-band 
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Fig. B.8. NGC 5047 V-band 



